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Abstract — Sampling rate is tlie bottleneck for spectrum sens- 
ing over multi-GHz bandwidth. Recent progress in compressed 
sensing (CS) initialized several sub-Nyquist rate approaches to 
overcome the problem. However, efforts to design CS reconstruc- 
tion algorithms for wideband spectrum sensing are very limited. 
It is possible to further reduce the sampling rate requirement and 
improve reconstruction performance via algorithms considering 
prior knowledge of cognitive radio spectrum usages. In this 
paper, we group the usages of cognitive radio spectrum into three 
categories and propose a modified orthogonal matching pursuit 
(OMP) algorithm for wideband spectrum sensing. Simulation 
results show that this modified OMP algorithm outperforms two 
modified basis pursuit de-noising (BPDN) algorithms in terms of 
reconstruction performance and computation time. 

I. Introduction 

Although the spectrum is almost fuUy allocated to various 
radio frequency (RF) technologies, it is not fully utilized. The 
spectrum measurement result in Fig. [I][l | shows that most of 
time, the utilization of the spectrum over 2.5 GHz bandwidth 
can be as low as 10%. This utilization is even lower in rural 
areas. Cognitive radio is the idea to increase the utilization. 
Spectrum sensing, the operation to find the unused spectrum 
holes in the whole spectrum, is essential to cognitive radio. 
Spectrum sensing for bandwidth over multi-GHz, however, 
posed a challenge to today's sampling technology. According 
to Nyquist sampling theorem, the sampling rate should be at 
least twice of the signal bandwidth in order to reconstruct the 
original signal. Current analog-to-digital converter (ADC) to 
reach multi-GHz sampling rate is hardly affordable. 

To reduce the sampling rate requirement, compressed sens- 
ing (CS) (HI, (]3] seems to be a natural approach. According 
to encouraging results from CS theory, the lowest sampling 
rate for perfect reconstruction is related to the sparsity, a.k.a. 
number of non-zeros, of the original signal. The sampling 
rate can therefore be reduced greatly compared with Nyquist 
rate. The CS theory has been successfully applied in many 
aspects of communications to reduce sampling rate [|4], Js], 
||6]. In cognitive radio, the spectrum is sparse in the frequency 
domain and CS has been applied to overcome the sampling 
rate problem. Several sampling structures have been proposed 
Q, |9|. These sampling structures have made CS based 
spectrum sensing possible, even in hardware prototypes lITOl . 
ifTTl . In this paper, however, we focus on reconstruction 
algorithms. 
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Fig. 1. Spectram utilization measured at Berkeley wireless research center. 



Standard CS algorithms have been used to reconstruct the 
original spectrum, such as basis pursuit (BP) liTSI . orthogonal 
matching pursuit (OMP) fT3\, etc. We propose to modify 
the standard OMP algorithm by using prior knowledge of 
cognitive radio spectrum usages to improve the reconstruction 
performance, thus further reducing the sampling rate require- 
ment. There are some cooperative wideband spectrum sensing 
algorithms modified upon standard CS algorithms llT4l . ifTSl . 
using shared knowledge of secondary users, under the assump- 
tion that there is channel to share information. In this paper, 
we do not assume secondary users to have channel for shared 
knowledge, and only prior knowledge for each individual 
secondary user is available. Very limited efforts have been 
made in modifying CS algorithms with prior knowledge for 
wideband spectrum sensing. An obvious prior knowledge of 
cognitive radio spectrum usages is fixed spectrum boundaries 
set by FCC for primary users. Such structural information 
can be used in CS algorithms as prior knowledge to improve 
the reconstruction performance. 'Mixed I2/I1 norm denoising 
operator' (MNDO) HH is an algorithm that modifies the BP 
problem formulation using such prior knowledge. Another BP 
based algorithm developed for magnetic resonance imaging 
(MRI), called modified BPDN (l3, has similar idea, and is 
introduced by us in wideband spectrum sensing for comparison 



TABLE I 

Three Categories of Cognitive Radio Spectrum Usage 

1. Spectram bands with fixed boundaries which are always accessed by 
primary users. Example: Local radio station signals, local TV signals, 
etc. 

2. Spectrum bands with fixed boundaries which are accessed by 
primary users not very often. Example: TV spectram available for 
white space device access 119]. 

3. Spectram bands with fixed boundaries which are partially and 
randomly accessed. Example: Cellphone signals. 



purpose. 

In this paper, we re-examine the prior knowledge of cog- 
nitive radio spectrum usages and propose a modified OMP 
algorithm. Similar to fTF|, the spectrum usages in a cognitive 
radio network can be grouped into three categories as listed 
in Table H] The modified OMP algorithm considers all three 
categories. We will show by simulation that the modified 
OMP algorithm has much better reconstruction performance 
in mean square error (MSE) than MNDO and modified BPDN, 
given the same amount of spectrum utilization. Moreover, our 
algorithm takes much less computation time compared with 
MNDO and modified BPDN. Such improvements are benefited 
from proper use of prior knowledge and the greedy nature of 
OMP algorithm. 

This paper is organized as follows. Section states the 
wideband spectrum sensing model under the CS framework. 
Section |III] introduces the MNDO, modified BPDN and our 
proposed modified OMP algorithm. Simulation results are 
shown in Section |IV] and conclusions are made in Section 

El 

II. Wideband Spectrum Sensing Model via 
Compressed Sensing 

A. Discrete Spectrum Model 

Suppose that the spectrum available for cognitive radio 
network has total bandwidth of B Hz. It is convenient to 
describe the spectrum in discrete domain using an x 1 
frequency domain sample vector f: 

{=[ h /2 ... fN f (1) 

where T is transpose, {/;} are signal values in frequency 
domain uniformly sampled over B with spacing B/N, and 
indices {i} are related to frequency locations. In noiseless 
case, if \fi\^ ^ 0, then the spectrum is occupied at the 
i-th frequency location. Otherwise, i represents the unused 
frequency location that is accessible to secondary users. The 
frequency sample vector f can be projected onto Nyquist rate 
time domain sample vector r via an inverse N x N Fourier 
projection matrix F^^: 

r = F-^i (2) 

According to the three categories mentioned in Table |I] 
fixed spectrum boundaries are known in prior Indices {bi} 
denote frequency boundaries and they can separate the whole 
spectrum into K consecutive subsections: 



ui = {l,2,...,6i} 

M2 = {&1 + 1,^1 + 2, ...,62} ^^-j 

UK = {bK-1 + 1, bK-1 + 2, N} 

Now, we define three category sets {5„} according to the 
following conditions: 

Sn ~ {ui \i G Category n} ,n — 1, 2, 3 

O = U 5n {S^ n Sj = 0, for i ^ j) (4) 

n 

where O is the universal set. 

According to the measurement results of spectrum utiliza- 
tion, we assume the following relation is valid 

\i\/N < 10% 
where |f | is the number of non-zeros in f. 

B. Compressed Sensing Based Spectrum Sensing Approach 

In the framework of CS ||20l . consider the sampling 
of an X 1 signal vector x = ^s, where ^ is an N x N 
dictionary matrix and s is an x 1 vector with L << N 
non-zero entries s, . It has been shown that x can be recovered 
with AI samples by projecting x via an M x N measurement 
matrix where L < M < N, and $ is incoherent with "S*. 
The M X 1 measurement vector y and the projection process 
can be represented as 

y = *X = ^'i'S (5) 

s can be reconstructed by BP yjl, which solves the 
following li norm convex optimization problem: 

i = arffmin ||s||i 

s " "1 (6) 

s.t. $\l/s = y 
where Ip norm of s for p > 1 is defined as: 

iisiip = (Ei^^i')' 

The signal can also be reconstructed by greedy algorithm 
like OMP 1 131, which will be described in detail. 

If we compare the system model in Equation (|2]l and 
Equation (|5]l, received time domain sample vector r can be 
considered as the signal vector x, inverse Fourier projection 
matrix F^^ can be considered as the dictionary matrix \I/ 
and frequency domain sample vector f can be considered 
as the sparse signal vector s. If a proper implementable 
measurement matrix 4> is designed to be incoherent with F^^, 
then the spectrum sensing problem can be modeled as the 
CS reconstruction problem, and sub-Nyquist rate sampling 
rate reconstruction can be possible via CS algorithms. As 
mentioned in Section |l] several efforts have been made for 
implementable measurement approaches in real world. In 
this paper, we focus on the reconstruction algorithm that 
can be applied to any of the sub-Nyquist rate measurement 
approaches. To emphasize the performance of reconstruction 



algorithms, we simply use an ill x Gaussian random matrix 
# to represent the measurement process together with inverse 
Fourier projection matrix F^^. The use of Gaussian random 
matrix can guarantee reconstruction performance (ISTI . The 
simplified spectrum sensing model is: 

y = *f (7) 

As a result, f can be reconstructed by solving the problem 

f = argmin llf |L 



s.t. €>f = y 



(8) 



For a more general model with additive white Gaussian 
noise (AWGN), we have 



y = $f + w 



(9) 



where w is M x 1 noise vector with normal distribution. The 
CS problem in ([Sj can be modeled using BPDN l,,12J : 

f = argmmi||y - *f 112+7 ||f 111 (10) 

where 7 is determined by noise level. (O is a special case of 
([Tol l by setting 7 to 0. 

In the following sections, we will compare different recon- 
struction algorithm performance using the simplified model in 
Equation (|9]l. 

III. Modified Orthogonal Matching Pursuit 

A. BP Based Reconstruction Algorithms Considering Prior 
Knowledge of Cognitive Radio Spectrum 

In order to consider the prior knowledge described in Table 
H we need to re-formulate the li norm problem in Equation 

Coll. 

MNDO lfT6l considered spectrum allocation boundaries as 
prior knowledge. The frequency domain sample vector f is 
divided into 'blocks' according to boundary information {hi}: 



f = 



/ij /hi J /fcl+lj /fc2 ' 



fl f2 

/fc/sT-l+ll ■••7 fbKl 



(11) 



MNDO re-formulates BP to a new convex optimization 
problem with boundary information: 



f = argmm(||fi||2 
s.t. ||y-*f 



< rj 



(12) 



where 77 is determined according to the noise variance. In 
this problem formulation, all values in the same block are 
correlated. They are prone to be zeros or non-zeros simul- 
taneously. All blocks are treated evenly, without considering 
the categories in Table U As will be shown in simulation in 
Section IIVI noisy signal will be reconstructed over multiple 
blocks even in noise-free scenario. 



Another reconstruction algorithm 'modified BPDN' has 
similar idea, and we propose it for cognitive radio spectrum 
sensing for comparison purpose. It groups the indices into two 
sets: 

• 'Dense' index set T — {i} where non-zeros are likely to 
exist. 

• 'Sparse' index set where non-zeros are not likely to 
exist and T"^ denotes the complement set of T. 

The modified BPDN algorithm models a convex optimiza- 
tion problem using the above prior knowledge is as follows. 



argmm- ||y 



*f|| 



- 7 l|fT<= 



(13) 



where 7 is determined by noise variance, denotes a sub- 
vector containing the elements of f with indices in T"^. In the 
solution of this convex problem, index sets in T are prone to 
have non-zero values, and index sets in T"^ are prone to have 
zeros values. 

We can use modified BPDN for wideband spectrum sensing. 
From Table |I] indices in Si are prone to have non-zeros 
while indices in ^2 and S3 are prone to have non-zeros. 
Therefore, it is straightforward to put 5*1 into dense index 
set T, while 5*2 and 6*3 into sparse index set T'^. However, 
if S2 and/or 5*3 have jointed subsets, some boundary indices 
will be invalid, resulting in invalid prior knowledge. Here 
we give a simple example of such invalid prior knowledge. 
Assume we have u„ = {bn-i + 1, 6,1-1 + 2...,6„} C S2 or 

S3, Urn = + l,&m-l + 2...,fe„i} C 5*2 Or S'3, with 

four boundary indices 6„_i, 5„, 6„i_i and If 6„ = 6,„_i, 
then there should be three valid boundary indices 5„_i, 6„ 
and bn+i- However, in the modified BPDN model, u„ Uum = 
{hn-i + 1, 6„, hn + 1, hn+i] C T"^ , and boundary index 
bn becomes invalid. As a result, the modified BPDN will 
reconstruct signal without considering the boundary index 5„ 
separating m„ and Um- 

B. Modifying the OMP Algorithm 

We have introduced two BP based reconstruction algorithms 
considering prior knowledge and their limitations. To over- 
come these, we propose a modified OMP algorithm, which 
will be able to deal with the three categories separately. 

OMP is a greedy algorithm. It reconstructs original signal 
by iterative search for non-zero indices and performs least- 
squares estimation of the values on the non-zero indices. The 
original OMP algorithm |13| is as follows: 

The iteration operation gives us the freedom to consider 
three categories separately. Since there must be non-zeros in 
^1, we can initialize the original Aq with 5*1. This modification 
benefits the reconstruction accuracy because occupied indices 
will always be counted. Then, during each iteration, if index 
At is found and if A* G Uj C 5*2 as well, all elements in Ui will 
be added to At. This modification benefits the reconstruction 
accuracy because only selected {ui} are counted. If index At £ 
Ui C 5*3, only At will be added, similar to the original OMP 
process. As a result, all three categories are treated separately. 
The modified OMP algorithm is summarized as: 



TABLE II 
The Original OMP Algorithm. 

Input: 

• An M X N matrix = * 

• An M X 1 sample vector y 

• Maximum number of iterations m 

• Error tolerance r] 
Output: 

• An estimate N X I vector f for the ideal signal 

• An index set At containing t elements from {1, 

• An M X I residual vector rest 

Procedure: 

1) Initialize the residual resg = y, the index set Aq = 
iteration counter t = 1. 

2) Find the index At that satisfies the following equation 



,N} 



At 



: arg max I (rest 
j = l,...,]V 



and the 



(14) 



3) 



4) 



where Bj denotes the j-th column vector of &, (a, b) denotes 
the inner product of two vectors a and b. If the maximum 
occurs for multiple indices, break the tie deterministically. 
Augment the index set At = At— i lj{-^t} and the matrix of 
chosen atoms &t = ©At> where ©At is a sub-matrix of 
containing columns with indices in At. 

Solve the least-squares problem to obtain a new signal estimate: 



xt = argmin |10tx - yjlj 

X 

5) Calculate the new residual: 



rest 



0txt 



(15) 



(16) 



6) If t < m or ||rest||2 > V' increment t, and return to Step 2. 

7) Xt is the estimated signal f , with non-zero indices at compo- 
nents Hsted in At. 



TABLE III 
The Modified OMP Algorithm. 

Input: 

• An M X A'^ matrix = * 

• An M X 1 sample vector y 

• Boundary infoiTnation {bi , 62 
{Su 82,83} 

• Maximum number of iterations m 

• Error tolerance r] 
Output: 

• An estimate N X 1 vector f for the ideal signal 

• An index set At containing t elements from {1, 

• An Af X 1 residual vector rest 
Procedure: 



bx}, {ui,U2, ■■■,uk} and 



.,N} 



y, the index set Aq = Si, the 
0Si and the iteration counter 



1) InitiaHze the residual reso 
matrix of chosen atoms 0i 
t = 1. 

2) Solve the least-squares problem in Equation il5\ to obtain a 
new signal estimate. 

3) Calculate the new residual using Equation | |16) . 

4) Increment t. 

5) Find the index At that satisfies Equation (14). 

6) Augment the index set At = At_i lJ{At}. If At G Ui C 52, 
let At = Atljui. 

1) Set the matrix of chosen atoms 0t = 0At ■ 

8) Solve the least-squares problem in Equation l llSt to obtain a 
new signal estimate. 

9) Calculate the new residual using Equation l ll6t . 

10) Return to Step 4 if t < m or ||rest!|2 ^ 

11) Xt is the estimated signal f, with non-zero indices at compo- 
nents listed in At. 



IV. Simulation Results 

In the simulation, we use a spectrum from 1 MHz to 
1000 MHz with approximately 10% utilization to test the 
performances of the reconstruction algorithms. There are 1000 
discrete frequency locations with 1 MHz resolution. Number 
of original siganl samples is = 1000. Assume boundary in- 
formation and category information as prior knowledge. Using 
Equation (O and (|4]l, we define three categories 81,82, 83 for 
demonstration purpose: 

n = {1, 2, 1000} = s*! u s'2 u 5*3 

81 = {ki,k2,...,h} 

where 



ki 


= {11,12,..., 


20} 


k2 


= {265,266, 


..,274} 


k3 


= {276,277, 


..,285} 


ki 


= {601,602, 


..,610} 




= {701,702, 


..,710} 



'S'2 — {^1, I2, ll5, hs, I 



19; 



, ^4o} 



where 

k ^ { 100 + 11 {i-l), 100 + 11 (i - 1) + 1, 100 + Hi - 

i = 1,2, ...,40 

^3 = 0/{5iU52} 

Note that this setup is just one example. It can be changed 
according to specific frequency ranges and frequency alloca- 
tions. 

When building up the spectrum utilization, all frequency 
locations in 5*1 are occupied, 10% frequency locations in ^2 
are randomly occupied, and 2% frequency locations in ^3 are 
randomly occupied. The total number of non-zeros is controled 
around 100. Fig. |2] illustrates one realization of the spectrum. 
We set the non-zero values to 1 for better illustration. The 
reconstruction algorithms can deal with arbitrary values. 
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Fig. 2. One realization of the spectrum utilization. 

We first demonstrate the sub-Nyquist rate reconstruction 
performance using Fig. |2] Let number of measurements 
M = 25%N = 250 and let * be 250 x 1000 Gaussian 



random matrix. Fig. |3] shows the results from modified OMP, 
MNDO and modified BPDN, under the noiseless case. We use 
normalized MSB as a metric to evaluate their performances: 



MSE : 



TABLE IV 

Averaged MSE for Noiseless Spectrum Realizations 





M = 250 


M = 500 


MNDO 


1.9975 


0.0810 


Modified BPDN 


0.1438 


7.920 X 10"^" 


Modified OMP 




7 X 10"-'^ 



MNDO has MSE « 1.8908. From Fig.|3] we can see that noisy 
spikes are generated over multiple frequency blocks within 
50 - 200 MHz and 1100 - 2000 MHz. Modified BPDN has 
MSE « 0.3331. Since boundary indices of jointed subsets 
become invalid, non-zero spikes can be found throughout the 
spectrum. Modified OMP, however, is able to reconstruct the 
original signal almost perfectly with MSE « 10" 



-30 



best performance among all of them, with MSE w 0.0386. 
The non-zero positions are perfectly found and errors are 
made from the fluctuation of the non-zero frequency values. 
Other algorithms have similar noisy signal reconstructed as the 
noiseless example. MNDO has MSE « 2.0007 and modified 
BPDN has MSE « 0.4248. 
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(c) Recovered by IVIodified BPDN 
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Fig. 3. (a) Recovered spectrum by modified OMP. (b) Recovered spectrum 
by MNDO. (c) Recovered spectrum by modified BPDN 

Now we use Monte Carlo simulation to compare the recon- 
struction performance for noiseless signal with M = 250 and 
M = 500. 500 different spectrum realizations are generated 
for the simulations. The averaged MSE performances for 
modified OMP, MNDO and Modified BPDN are in Table [IVl 
It can be seen that when M — 250, only modified OMP 
can provide almost perfect reconstruction with MSE w 1'^°. 
When M — 500, modified BPDN can reach almost perfect 
reconstruction with MSE w 10"^", while MNDO still has a 
few eiTors with MSE ?» 0.081. 

Next we show the example of the algorithms' reconstruction 
ability under AWGN. The original noisy spectrum is shown 
in Fig. m with signal-to-noise ratio (SNR) set to 15 dB. 
With M — 250, the reconstructed spectrum using these 
algorithms are shown in Fig. |5] Obviously, modified OMP has 
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Fig. 4. One realization of the noisy spectrum utilization. SNR = 15 dB. 

Monte Carlo simulation is performed in the noisy signal re- 
construction. We use M = 250 and M = 500 for simulations. 
SNR is changed from 5 dB to 20 dB with 5 dB increments. Fig. 
|6]and Fig. |2]show the simulation results. For each SNR plot, 
500 noisy signal realizations are used. Obviously, modified 
OMP outperforms other two BP based algorithms in AWGN 
as well. Better MSE performance for all algorithms can be 
achieved by increasing SNR or M. 

The final test of the reconstruction algorithms is the com- 
putation time. All algorithms are written in Matlab and run 
on a computer with 2.8 GHz Pentium 4 CPU. The code for 
modified OMP is written in basic Matlab functions. The codes 
for MNDO in (O and modified BPDN in ([T3]l are written 
based on cvx, a package for specifying and solving convex 
programs Il22l . ||23i . To get the results in Fig. |3] modified OMP 
takes about 0.05 seconds, MNDO takes about 12.83 seconds, 
and modified BPDN takes about 5.89 seconds. Obviously, 
modified OMP is much faster 

V. Conclusion 

By grouping cognitive radio spectrum usages into three 
categories, we propose a modified OMP for cognitive radio 
wideband spectrum sensing. For comparison purpose, we 
introduce modified BPDN for wideband spectrum sensing 
as well. The modified OMP is compared with MNDO and 
modified BPDN, which are BP based. We show by simulation 
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Fig. 5. (a) Recovered spectrum by modified OMP. (b) Recovered spectrum 
by MNDO. (c) Recovered spectrum by modified BPDN. 




Fig. 6. M = 250, MSE performances for all algorithms under AWGN. 



and conclude that, by grouping cognitive radio spectrum 
usages into three categories, the proposed modified OMP 
outperforms other two BP based algorithms in both accuracy 
and computation efficiency. The improvements in accuracy and 
speed lie in the proper use of prior knowledge and the greedy 
aspects of the OMP. 

Based on the OMP framework, more work is encouraged 
to be done to explore new algorithms in cooperative spectrum 
sensing, using reasonable prior knowledge to further improve 
the performance of CS based wideband spectrum sensing. 
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Fig. 7. M = 500, MSE performances for all algorithms under AWGN. 



Acknowledgment 

This work is funded by National Science Foundation 
through grants (ECCS-0901420), (ECCS-0821658) and Office 
of Naval Research through two contracts (NOOO 14-07 -1-0529, 
N00014-1 1-1-0006). 



References 



[1] 



S. Mishra, D. Cabric, C. Chang, D. Willkomm, B. van Schewick, 
S. Wolisz, and B. Brodersen, "A real time cognitive radio testbed 
for physical and link layer experiments," in First IEEE International 
Symposium on New Frontiers in Dynamic Spectrum Access Networks, 
pp. 562 -567, 8-11 2005. 
[2] D. Donoho, "Compressed sensing," IEEE Transactions on Information 

Tlieory, vol. 52, no. 4, pp. 1289-1306, 2006. 
[3] E. Candes, "Compressive sampHng," in Proceedings oftlte International 

Congress of Matliematicians, vol. 3, pp. 1433-1452, 2006. 
[4] P. Zhang, Z. Hu, R. Qiu, and B. Sadler, "A compressed sensing based 
ultra-wideband conmiunication system," in Proc. IEEE International 
Conference on Communications ICC '09, pp. 1-5, 2009. 
[5] Z. Wang, G. R. Arce, J. L. Paredes, and B. M. Sadler, "Compressed 
detection for ultra-wideband impulse radio," in IEEE Sig. Proc. Adv. 
Wireless Comm. (SPAWC), 2007. 
[6] J. Paredes, G. R. Arce, and Z. Wang, "Ultra-wideband compressed 
sensing: Channel estimation," IEEE J. Select. Topics Signal Proc, vol. 1, 
pp. 383-395, Oct. 2007. 
[7] Z. Tian and G. Giannakis, "Compressed sensing for wideband cognitive 
radios," in IEEE International Conference on Acoustics, Speecli and 
Signal Processing, vol. 4, pp. 1357-1360, 2007. 
[8] Y. Polo, Y. Wang, A. Pandharipande, and G. Leus, "Compressive wide- 
band spectrum sensing," in IEEE International Conference on Acoustics, 
Speech and Signal Processing, pp. 2337 -2340, 19-24 2009. 
[9] M. Mishali and Y. Eldar, "From theory to practice: Sub-nyquist sampling 
of sparse wideband analog signals," IEEE Journal of Selected Topics on 
Signal Process, vol. 4, pp. 375-391, Apr 2009. 
[10] T. Ragheb, J. Laska, H. Nejati, S. Kirolos, R. Baraniuk, and Y. Massoud, 
"A prototype hardware for random demodulation based compressive 
analog-to-digital conversion," in 5Lst Midwest Symposium on Circuits 
and Systems, 2008, pp. 37^0, IEEE, 2008. 
[11] M. Mishali, Y. Eldar, O. Dounaevsky, and E. Shoshan, "XampHng: 
Analog to digital at sub-nyquist rates," lET Journal of Circuits, Devices 
and Systems, 2009. to appear. 
[12] S. Chen, D. Donoho, and M. Saunders, "Atomic decomposition by basis 

pursuit," SIAM review, vol. 43, no. 1, pp. 129-159, 2001. 
[13] J. Tropp and A. Gilbert, "Signal recoveiy from random measurements 
via orthogonal matching pursuit," IEEE Transactions on Information 
Tlieory, vol. 53, no. 12, pp. 4655^666, 2007. 



[14] Z. Zhang, H. Li, D. Yang, and C. Pei, "Space-time bayesian compressed 
spectrum sensing for wideband cognitive radio networks," in 2010 IEEE 
Symposium on New Frontiers in Dynamic Spectrum, pp. 1-11, 6-9 2010. 

[15] Y. Wang, A. Pandharipande, Y. Polo, and G. Leus, "Distributed com- 
pressive wide-band spectrum sensing," in Information Theory and Ap- 
plications Workshop, pp. 178 -183, 8-13 2009. 

[16] Y. Liu and Q. Wan, "Compressive wideband spectrum sensing for fixed 
frequency spectrum allocation." Arxiv preprint arXiv:1005. 1804, 2010. 

[17] N. Vaswani and W. Lu, "Modified basis pursuit denoising (modified- 
bpdn) for noisy compressive sensing with partially known support," 
in IEEE International Conference on Acoustics, Speech and Signal 
Processing, 2010. 

[18] S. Haykin, "Cognitive radio: brain-empowered wireless communica- 
tions," IEEE Journal on selected areas in communications, vol. 23, no. 2, 
pp. 201-220, 2005. 

[19] S. Jones, T. Phillips, H. Tuyl, and R. Weller, "Evaluation of the 
performance of prototype TV-Band white space devices, phase n," tech. 
rep., FCC OET Report 08-TR-1005, Oct. 15 2008. 

[20] E. Candes and T. Tao, "Near-optimal signal recovery from random 
projections: Universal encoding strategies?," IEEE Transactions on 
Information Theory, vol. 52, no. 12, pp. 5406-5425, 2006. 

[21] E. Candes, J. Romberg, and T. Tao, "Robust uncertainty principles: Exact 
signal reconstruction from highly incomplete frequency information," 
IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489-509, 
2006. 

[22] M. Grant and S. Boyd, "CVX: Matlab software for disciplined convex 
programming, version 1.21." http://cvxrcom/cvx, Jul. 2010. 

[23] M. Grant and S. Boyd, "Graph implementations for nonsmooth convex 
programs," in Recent Advances in Learning and Control (V. Blon- 
del, S. Boyd, and H. Kimura, eds.). Lecture Notes in Control and 
Information Sciences, pp. 95-110, Springer- Verlag Limited, 2008. 
http://stanford.edu/~boyd/graph_dcp.html. 



